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ABSTRACT 



Using the Moon to occult the Sun, the Clementine spacecraft used its navigation cameras to 
map the inner zodiacal light at optical wavelengths over elongations of 3 < e < 30° from the Sun. 
This surface brightness map is then used to infer the spatial distribution of interplanetary dust 
over heliocentric distances of about 10 solar radii to the orbit of Venus. The averaged ecliptic 
surface brightness of the zodiacal light falls off as Z(e) oc e - 2 - 45±0 05 which suggests that the dust 
cross-sectional density nominally falls off as a(r) oc r -!- 4 5±0-05_ interplanetary dust also has 
an albedo of a ~ 0.1 that is uncertain by a factor of ~ 2. Asymmetries of ~ 10% are seen in 
directions east-west and north-south of the Sun, and these may be due the giant planets' secular 
gravitational perturbations. 

We apply a simple model that attributes the zodiacal light as due to three dust populations 
having distinct inclination distributions, namely, dust from asteroids and Jupiter-family comets 
( JFCs) having characteristic inclinations of i ~ 7° , dust from Halley-type comets having i ~ 33° , 
and an isotropic cloud of dust from Oort Cloud comets. The best-fitting scenario indicates that 
asteroids + JFCs are the source of about 45% of the optical dust cross-section seen in the ecliptic 
at 1 AU, but that at least 89% of the dust cross-section enclosed by a 1 AU radius sphere is of a 
cometary origin. Each population's radial density variations can also deviate somewhat from the 
nominal a(r) oc r -1 ' 45 . When these results are extrapolated out to the asteroid belt, we find an 
upper limit on the mass of the light-reflecting asteroidal dust that is equivalent to a 12 km asteroid, 
and a similar extrapolation of the isotropic dust cloud out to Oort Cloud distances yields a mass 
equivalent to a 30 km comet, although the latter mass is uncertain by orders of magnitude. 
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1. Introduction 

Interplanetary dust is of considerable interest since these grains represent samples of small 
bodies that formed in remote niches throughout the solar system. Dust grains are liberated when 
the rocky asteroids collide and when icy comets sublimate during a close approach to the Sun, 
and this dust is transported throughout the solar system by solar radiation forces. If it is granted 
that the information carried by this dust is indeed decipherable, then samples of this dust tell us 
about the conditions in various parts of the solar nebula from which asteroids, comets, and planets 
subsequently formed. In particular, asteroidal dust tells us of nebula conditions at the boundary 
between the terrestrial and the giant-planet zones in the solar nebula. And since the long-period 
comets from the Oort Cloud formed amongst the giant planets, the mineralogy of their dust is 
indicative of conditions over a vast swath of the solar nebula between ~ 5-30 AU. Information 
about the outer reaches of the solar nebula is also carried by dust generated by the shorter-period 
Jupiter-family comets that likely formed in the Kuiper Belt beyond ~ 30 AU. 

This dust is also of dynamical interest since the spatial density of these 'trace particles' allows 
one to simultaneously assess the relative strengths of Poynting-Robertson drag (which drives dust 
sunwards), the planets' gravitational perturbations (which disturbs dust orbits and sometime con- 
fines dust at resonances), and mutual collisions (which fragments and destroys dust). However a 
deeper understanding of this dust first requires knowledge of the abundance and spatial distribution 
of asteroidal and cometary dust grains, both of which are the subject of this investigation. 

Whipple (1955) performed one of the earliest assessments of the various sources of interplane- 
tary dust. This analysis is essentially a mass-budget that compares the rates at which comets and 
asteroids produce dust to the rates at which collisions and Poynting-Robertson drag destroy dust. 
From the very limited data on comets, asteroids, and interplanetary dust that were available at 
the time, it was concluded that at least 90% of interplanetary dust is of cometary origin (Whipple 
1955, 1967). This view prevailed for the next three decades until the Infrared Astronomical Satellite 
(IRAS) discovered the asteroidal dust bands (Neugebauer et al. 1984; Dermott et al. 1984). These 
dust bands are clearly produced by asteroid families, so these observations demonstrate that as- 
teroids are also significant contributors to the interplanetary dust complex, with models indicating 
that ~ 30-40% of the outer zodiacal light is due to asteroidal dust (Dermott et al. 1994; Liou et 
al. 1995). Asteroidal and cometary dust have also been collected in the Earth's stratosphere by U2 
aircraft (Brownlee et al. 1993). Atmospheric entry velocities can be inferred from these grains, and 
about 80% of the dust in this sample have low entry velocities consistent with asteroidal orbits. 
However this finding should be regarded as an upper limit on the true abundance of asteroidal dust 
in the ecliptic since (i.) dust released from low-inclination Jupiter-family comets (JFCs) can also 
have low entry velocities (see Section 4.2.2), and (ii.) the Earth's gravitational focusing naturally 
selects for low-velocity dust from asteroids (and JFCs as well) over dust from the higher-inclination 
Halley-type and Oort Cloud comets. 

In order to assess the abundances and spatial distributions of asteroidal and cometary dust in 
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the inner solar system, the following analyzes images of the zodiacal light that were acquired by 
the Clementine spacecraft while in lunar orbit. Clementine orbited the Moon for about two months 
in early 1994, ostensibly to study the lunar surface. But a secondary objective of this mission was 
to image the inner zodiacal light using Clementine's wide-angle navigation cameras. While the 
Sun was in eclipse behind the Moon, the Clementine star tracker cameras acquired hundreds of 
images of the zodiacal light over elongations that span the orbit of Venus down to about 10 solar 
radii. As this is the first scientific application of a star tracker camera, the instrument and its 
optics are described in detail in Section 2 and Appendix A. Section 3 and Appendix B describe 
the observations and data reductions, also in some detail due to several artifacts present in the 
data. However the reader uninterested in these particulars can skip directly to Section 4 where the 
interplanetary dust model is described and applied. Results are then summarized in Section 5. 



2. The Star Tracker Camera 

The zodiacal light images studied here were acquired by the Clementine spacecraft's two star 
tracker cameras. A star tracker is a simple, light-weight, low-power camera designed to acquire 
wide-angle CCD images of star fields. The spacecraft's two star tracker cameras are designated 
A and B, and nearly all of the data examined here were acquired by star tracker B. The principle 
purpose of the star tracker is to aid spacecraft navigation; by comparing the observed star fields to 
an onboard star atlas, the spacecraft can continuously monitor its orientation. It should be noted 
that a high-quality photometric imaging capability was not a design criterion for this camera. 
Nonetheless, our close inspection of the data shows that this camera can be used to obtain high- 
fidelity images once a number of instrumental artifacts are removed from the data. These data 
reductions are described in detail in Section 3. 

A simplified schematic of the star tracker optics is shown in Fig. 1, and a more detailed 
description of the instrumentation may be found in Lewis et al. (1991) and Kordas et al. (1995). 
The camera's principle components are a spherical lens, a fiber optic, and a CCD detector. The 
focal point of this lens is at its center. As Fig. 1 shows, incident light entering the lens from the 
left forms an image at the opposite side of the lens, and the fiber optic pipes this light to the CCD. 
Appendix A shows how to map the CCD's (x, y) coordinates for every pixel into equatorial and 
ecliptic coordinates. 

The CCD detector is a Thomson TH7883 array of 384 x 576 pixels. The camera's angular 
field of view is 29.3° x 44.6° and the plate-scale at the optical axis is 0.0756 degrees /pixel. This 
camera's point spread function has a full width at half maximum of about 2 pixels ~ 0.15°. No 
filter was used in this camera. The detector has a peak quantum efficiency of about 45% at a 
wavelength of A ~ 8000 A (Lewis et al. 1991). Figure 2 gives the camera's relative instrumental 
response and shows that the camera is sensitive to wavelengths of 5000 < A < 9000 A. Figure 2 
also shows the equivalent square bandpass (e.g., one having the same area under the curve as the 
observed instrumental response) that has a spectral width of AA = 3490 A and a mean wavelength 
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A=7370 A. 

It should be noted that this camera also has two serious handicaps. The first is that the 8-bit 
CCD has a dynamic range of only 256. However the star tracker images were acquired at exposures 
that differ by up to a factor of 14, so the total dynamical range of the zodiacal map produced here 
is 14 x 256 ~ 3600. Another problem is the absence of a shutter in the camera which results in the 
CCD being continuously exposed as the array is read out along the detector's columns. This tends 
to redistribute the flux from all sources along the CCD columns. However as Section 3 shows, this 
effect is reversible, and 'destreaked' data may be recovered from the raw data itself. 

3. Observations and Data Reductions 

The Clementine spacecraft was in an elliptical polar orbit about the Moon from February 22 
through May 4 of 1994, after which it left the Moon for an encounter with the near-Earth asteroid 
Geographos and was subsequently lost due to a software failure. But during the final six weeks in 
orbit about the Moon, the star tracker cameras repeatedly imaged the inner zodiacal light while the 
Moon occulted the Sun. Numerous images were acquired during an orbit about the Moon, either 
just prior to sunrise or just after sunset. Each batch of images are identified by an orbit number 
that is simply the number of lunar orbits that Clementine had achieved to date. Because these 
observations were acquired during a six week interval, the longitude of the camera's line-of-sight 
to the Sun changed considerably due to the heliocentric motion of the Earth-Moon system. Table 
I lists orbit numbers, observation dates, each observation's heliocentric ecliptic longitude, and total 
exposure times for the subset of the data that are examined here. The camera's lines of sight 
through the ecliptic during different orbits are also shown in Fig. 3. 

All of the raw Clementine data examined here are archived at the National Space Science 
Data Center and may be obtained at the URL 

http://nssdc.gsfc.nasa.gov/planetary/lunar/clementinel.html, and the flatfield used to process the 
star tracker images is available from the authors. 

The star tracker camera usually acquired a sequence of about 40 or so images during each orbit 
of the spacecraft. The exposure times for every image acquired during a sequence usually cycled 
between 0.05, 0.1, 0.2, 0.4, and 0.7 seconds. By cycling the exposure times in this manner, the 
camera's dynamic range was increased by a factor of 14 and difficulties due to image-saturation in 
the brighter parts of the zodiacal light were mitigated. A typical raw image is shown in Fig. 4A, 
which is a 0.4 sec exposure acquired during orbit 193. Although the Sun is well behind the lunar 
limb, the Moon is partly illuminated by sunlight reflected by an Earth that is outside of the field 
of view. The bright object left of the Moon is a saturated Venus. In fact, Venus is so bright that 
the signal accumulating in the pixel at Venus has bled into the nearest 10 or so pixels. 

Note also the bright streak at Venus in Fig. 4A as well as a broader but dimmer streak running 
through the core of the zodiacal light. These streaks are a consequence of reading the CCD array in 
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a shutterless camera. The camera electronics reads the CCD array by shifting the contents of every 
pixel down along the CCD's columns that, in Fig. 4, run left -to-right. As a row of "logical" pixels 
shifts off the bottom of the array, their values are recorded, zeroed, and then that row shifts back 
into the CCD's top row. However the CCD is always exposed during this process, so reading out 
the CCD array causes every logical pixel to receive additional signal from all parts of the sky that 
subtend that pixel's CCD column. This results in an image that appears to have streaks running 
along the columns (see Fig. 4A). Nonetheless, the time to transfer the contents of one pixel to the 
next, 94.4 /xsec (Kordas et al. 1995), is a fixed quantity, so the intensity of each column's streak can 
be inferred and removed from the raw data itself using the destreaking algorithm given in Zook et 
al. (1997). However this algorithm fails whenever a CCD column contains one or more saturated 
pixels. In this case, the pixels along the entire column are flagged and are subsequently ignored 
for the remainder of the analysis. Pixels polluted by Venus also disturb the destreaking algorithm, 
so they are first replaced with a local average of the zodiacal light prior to destreaking and are 
subsequently discarded (see Fig. 4B). 

In principle the dark current should be subtracted before destreaking an image. However the 
lack of a shutter makes it impossible to directly measure the dark current from these streaked 
images. Nonetheless, images acquired during orbit 66 show a very dark lunar surface that is 
shadowed from both the Sun and the Earth, so the flux observed in lunar surface provides a first 
estimate of the dark current. An aperture is placed on the dark portion of the Moon and an 
initial estimate of the dark current / is obtained for every image in the sequence. This current 
is subtracted from each image which is then destreaked. The residual flux Sf in that aperture is 
then examined, and automated software then revises the estimated dark current / appropriately 
and this cycle repeats until 5f has relaxed to zero. A similar algorithm is also used to subtract 
the dark current from all of the other images acquired during different orbits. However these 
images generally show a lunar surface that is either slightly or wholly illuminated by earthshine, 
so the flux measured in the lunar aperture represents the dark current plus a nonzero offset Sf. In 
this case, the above algorithm iteratively subtracts the dark current / and destreaks each image 
until a predetermined residual flux Sf > is achieved. The value for Sf appropriate for each 
image-sequence is determined later by comparing images that overlap the orbit 66 field; see below. 

The construction of the star tracker flatfield is described in Appendix B. Each image is 
flatfielded and pixels that subtend the Moon are flagged and discarded. Next, small shifts to the 
images are applied as needed so that the stars seen in an image-sequence appear stationary. A 
single averaged "master" image is then formed from the image-sequence using only the good pixels 
that were not previously flagged as bad. Figure 4B shows the master image for orbit 193; the 
data-gaps correspond to the Moon as well as pixels polluted by Venus. 

Two additional faint artifacts become evident upon close inspection of this longer-exposure 
image. The first is that every eighth column (which runs left -right in Figure 4B) is slightly darker 
than its neighbors. It is only evident at the outer edges of the images where the zodiacal light is quite 
faint, and its effect is barely discernible in Fig. 4B. This is probably due to a slight inhomogeneity 
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in the dark-current across the CCD. However this faint striping is of little consequence since its 
magnitude is comparable to the noise in the data. Another faint dark stripe can also be seen in 
the rows that subtend Venus; its magnitude is roughly twice the pixel-to-pixel noise in the image. 
The cause of this stripe is unknown, and it is only seen in images containing a deeply saturated 
Venus. 

Figure 4C shows a textbook quality image of the zodiacal light which has the Moon and Venus 
pasted back in. 

The camera's plate-scale and its pointing are determined from the handful of bright field stars 
that are identified in each of the master images; see Appendix A for details. Using each star's 
observed (x, y) coordinates, their known equatorial coordinates (a, 5), and Eqns. (Al) and (A6), 
a plate-scale of p — 0.0756 degrees/pixel is obtained. The lens coordinates (9,<j>) for each star arc 
then computed (see Fig. 1 and Eqns. Al), and Eqns. (A3) are solved for (a ,5 ), which are the 
equatorial coordinates for the pixel at the camera's optical axis, and r, which is the angle between 
the CCD's y axis and equatorial north. With these quantities known, equations (Al) and (A3) can 
now be used to compute equatorial coordinates for every pixel in each master image, and Eqs. (A5) 
are used to rotate these coordinates into geocentric ecliptic longitude and latitude (A, 6), as well 
as the longitude of each pixel relative to the Sun, A — A . 

With the pointing for every master image known, it is now possible to determine the unknown 
offsets 5f for the remaining master images that have not yet had their dark current properly 
subtracted. By examining those fields that overlap the orbit 66 field, it is straightforward to 
estimate the small offsets 5f that yield a mutually consistent surface brightness in the overlapping 
regions. This process is then repeated for the remaining adjacent fields until all offsets for all 
images have been determined. With these new offsets in hand, the entire data-reduction cycle 
(dark current subtraction, destreaking, flatfielding, and offset determination) is repeated until no 
further changes in the offsets are required. 

The final step is to merge all of the master images into a single wide-angle mosaic of the 
inner zodiacal light. This results in the 60° x 60° mosaic shown in Fig. 5. This image is formed 
by mapping the intensity of every good pixel in all of the master images into the corresponding 
2x2 box of pixels in the mosaic, which smooths the mosaic over an angular scale of 2p = 0.15°. 
Also recall that these images were acquired over a six week interval, so some stars are seen more 
than once as they drift to the right with time due to the heliocentric motion of the spacecraft. 
Consequently, several planets are also seen at multiple longitudes: Saturn is barely discernible at 
about 11° west of the Sun, Mars appears at 16° and 18° west, Saturn again at 19°, and Mercury 
at 27° west. If the saturated Venus were not already clipped from these images, it would inhabit 
the data-gap at 20° east of the Sun, 

The observed intensity of (5 Hydri is used to calibrate these data. This G2IV star is the only 
bright object in these images having a solar-type spectrum. This star has B — V and U — B colors 
that are solar to within 0.05 magnitudes, has an apparent V magnitude of m+ = 2.80, and has 
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an instrumental intensity of 1+ = 770 ± 60 counts/sec. Note that this intensity is obtained from 
images that are flatfielded using the "point-source" flatfield that is described in Appendix B. One 
common brightness unit in zodiacal light observations is the mean solar brightness B & = 
which is the intensity of the Sun 7 Q divided by the solid angle of the Sun Q & = 0.223 deg 2 . Since 
h = J lO-°- 4 ( m *- m ©) for a solar-type star, it follows that 

1 coun t / S ec/ P ixel = 10^-—)^ (55^)" ^ (D 

where tuq = —26.78 is the apparent visible magnitude of the Sun and Q p = p 2 = 5.72 x 10~ 3 deg 2 
is the solid angle of a pixel at the optical axis. Another common unit is SW@ = 4.33 x lO -16 !?©, 
which is the intensity of a tenth magnitude solar type star distributed over a square degree. Thus 
1 count /sec/pixel = (7.5 ± 0.6) x 1O _14 -B = (170 ± 10)S10q. The seemingly large uncertainty of 
8% in this calibration constant is due to (a) (3 Hydri's short exposure time of only 0.6 sec, and (b) 
noise in the flatfield — see Appendix B. 

However the relative uncertainties in the mosaic image, Fig. 5, vary across of the field due 
to the different exposure times of the various master images (see Table I). The fields west of 
Sun, which were acquired during orbits 66 and 110, had very short exposure times, so the western 
side of the mosaic is considerably noisier than the eastern side. Uncertainties in the dark current 
subtraction are ~ 3 counts/sec ~ 2 x 1O~ 13 -B ~ 5OOS , 1O . However this uncertainty is significant 
only at the outer edges of Fig. 5 where it can be as much as 50% of the signal there. We also note 
that nearly all of the images acquired during orbit 164 had a bright Earth in its field of view, so 
these images have considerable amounts of scattered light in them. This is the field just north of 
the Sun in Fig. 5, and this polluted zone lies at elongations of <E> > 10° north and north-northwest 
of the Sun. This is the only field acquired by star tracker A for which the flatfield is unavailable. 
We have elected to process this field using the flatfield from star tracker B, and it is included in 
the mosaic Fig. 5 solely for the purpose of filling an otherwise large datagap. The light-polluted 
portion of this field is not used in the subsequent analysis. 

Faint, diffuse background light from the galaxy also contaminates Fig. 5. However this was 
minimized by observing at an epoch when the sunward lines of sight were at the highest possible 
galactic latitudes of 30° < j3 g < 90° (see Fig. 3). The surface brightness of the galaxy was measured 
by Pioneer 10 while at heliocentric distances r > 2.8 AU where the zodiacal contribution is negli- 
gible; at latitudes f3 g > 30° the galactic surface brightness is Z g < 90 S10q (Leinert et al. 1998) 
at the southern edge of Fig. 5, and it decreases to the north. However this flux is substantially 
smaller than the uncertainty in the dark current subtraction and is neglected here. 

The integrated intensity of the light seen in Fig. 5 is / = 4.8 x 10~ 8 /q which corresponds to a 
visual magnitude my = —8.5. This makes the zodiacal light the second brightest object in night 
sky, the first being the full Moon having my = —12.7 and the third being Venus with my = —4.6 
at its brightest. 



-10- 



4. A Simple Model of the Interplanetary Dust Complex 



In order to extract the gross properties of the observed dust seen in Fig. 5, a simple model 
that is quite common in the literature shall be fitted to these data [c.f., Leinert (1975)]. The 
model assumes that the dust density varies radially as a power-law with heliocentric distance and 
that the cloud is axially symmetric. This model also assumes that the center of the cloud is at 
the Sun and that its midplane is in the ecliptic. Although none of these assumptions are actually 
correct in detail, they are sufficiently good for our purposes. In this case the spatial density of dust 
cross-section a can be written as a function of the heliocentric distance r and heliocentric ecliptic 
latitude (3: 



where r± = 1 AU is a reference distance, o\ = a(r\,0) is the dust cross-section density in the 
ecliptic at r = r±, and h{(3) describes how the dust density falls off with ecliptic latitude. 

The surface brightness of the sunlight reflected by this dust distribution is proportional to 
a(r, (3) multiplied by a light -scattering function and integrated along an observer's line-of-sight. 
The flux density of sunlight that is reflected by dust in a small volume element dV is dF = 
a(r, [3)<S>(p)(L Q /4irr 2 )dV/ A 2 (Lester et al. 1979) where dV = ttA 2 dA and Q is the solid angle of 
the volume element as seen by an observer a distance A away; see Fig. 6 for the definition of all 
the geometric quantities used here. The scattering phase function <&(</?) is related to the phase law 
ip{p) via <£(</?) = (o,/ir sr)ifj(ip) where a is the dust geometric albedo and <p is the scattering angle. 
Note that this formulation is valid only in the geometric optics limit, which is indeed the case since 
the bulk of the dust cross-section is contributed by grains having sizes ~ 10-100 /j,m (Grim et al. 



The surface brightness of the zodiacal light is thus Z = J dF/Q integrated over < A < oo. 
Noting that rjr\ = sine/ sint/? where e is the elongation of the line of sight having a geocentric 
ecliptic latitude and longitude (9,<fi) relative to the Sun, then cose = cos 4> cos 6, A/r\ = sin(</? — 
e)/ sin(^, dA/r\ = sineiii^/ sin 2 ip, and so the surface brightness can be recast as an integral over 
the scattering angle ip (Aller et al 1967; Giese and Dziembowski 1969): 




(2) 



1985). 



Z{9,4>) = 




B Q f \l>(<p)h(J3(<p))sm''(<p)d<p 



J c 



(3) 



where Q Q = 0.223 deg 2 = 6.80 x 10" 5 sr is the solid angle of the Sun and B Q = L Q /4Trr 2 Q Q is the 
mean surface brightness of the Sun where L© is the solar luminosity. Note also that the (3 in Eq. 
3 depends on the scattering angle ip through sin (3 = sin((p — e) sin 0/ sine. 
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4.1. Radial variations 

For a line of sight in the ecliptic, /t(0) = 1 and Eq. (3) becomes a simple integral over the 
phase law ip(tp). Two very different empirical phase laws are shown in Figure 7. The upper curve 
was constructed by Lamy and Perrin (1986), and it exhibits a very strong forward scattering peak 
(i.e., ip diverges as <p — > 0) as might occur due to the diffraction of sunlight by dust larger than 
a wavelength. The lower phase law is from Hong (1985); although this law does not show any 
forward scattering, it does exhibit a mild backward scattering peak at <p ~ 180° as is required of 
any phase law in order to reproduce the gegenshein. Despite the very different forms for VKv 7 ); 
both phase laws are very able at reproducing a varied suite of other zodiacal light measurements 
that were acquired over a wide range of elongation angles e (Hong 1985; Lamy and Perrin 1986). 
Consequently, our results are remarkably insensitive to the choice of the phase law. Regardless of 
whether one adopts Hong's backscattering phase function or Lamy and Perrin's forward scattering 
phase law, a numerical integration of the Eq. (3) yields a surface brightness Z(e) that, over our 
observation interval 2 < e < 30°, is largely indistinct aside from a numerical factor of ~ 1.6. 
This insensitivity to the details of ip(<P) 1S due to the fact that the dominant contribution to the 
surface brightness integral is by dust in the vicinity of 92 ~ 90°, i.e., dust nearer the Sun along the 
line-of-sight. 

This particular behavior also means that the integral in Eq. (3) is quite insensitive to the lower 
integration limit for the range of elongations e < 30° that are considered here. In this case the line 
of sight integral evaluates to ~ 0.83 when the Hong phase law is used and ~ 1.3 when the Lamy 
and Perrin law is used. Eq. (3) then simplifies to 

Z( e )^(2.3±O.5)xlO- 5 ^^5 (4) 

sin T e 

where the error in the coefficient indicates the uncertainty in the phase law. We also note that 
about 90% of the light seen in the ecliptic at elongations e < 30° is contributed by dust orbiting 
interior to 0.6 AU. 

The radial power law v is now readily obtained from profiles of the zodiacal light's ecliptic 
surface brightness. East-west and north-south profiles are shown in Fig. 8, and a power-law fit to 
the averaged east-west profile yields 

Z(e) - (I-? ± 0-2) x 10-" ( ) 

Z i e ) ~ • 2.45±0.05 ; W 

so a comparison with Eq. (4) shows that v = 1.45 ± 0.05 and ao\r x = (7.4 ± 1.8) x 10~ 9 . The 
uncertainty in the former quantity includes the statistical variations of the data seen in Fig. 8 while 
the latter quantity also includes the 8% uncertainty in the calibration and the uncertainty due to 
the possible choices for the phase law. Note that the power law v reported here is slightly steeper 
than that inferred from the data obtained by the Helios spacecraft (Leinert et al. 1981) and the 
Cosmic Background Explorer (COBE) spacecraft (Kelsall et al. 1998), with the disagreements at 
the 2cr level. 
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Figure 8 also reveals an asymmetry in the surface brightness of the zodiacal light north/south 
of the ecliptic, as well as an asymmetry east /west of the Sun. Ratios of the north/south and 
east /west surface brightness profiles are plotted in Fig. 9 which shows that the northern hemisphere 
gets steadily brighter with elongation relative to the southern hemisphere, as does the eastern ansa 
relative to the western ansa. Asymmetries such as these have been attributed to the giant planets' 
secular gravitational perturbations which can organize the longitudes of the dust grains' perihelia 
and nodes (Wyatt et al. 1999). The north-south asymmetry seen here is likely the same asymmetry 
previously observed by the Helios 1 and 2 which detected a i = 3° tilt between the midplane of 
the inner zodiacal light and the ecliptic (Leinert et al. 1980). The node of this symmetry plane has 
a longitude of $7 = 87° which, as Fig. 3 shows, is largely perpendicular to the Clementine lines- 
of-sight, and this particular viewing geometry will make one hemisphere slightly brighter than the 
other. Gravitational perturbations by giant planets can also displace the zodiacal cloud's center of 
light radially away from the Sun, which results in a phenomenon known as pericenter glow (Wyatt 
et al. 1999); such perturbations may be responsible for the east-west asymmetry seen in Figs. 8-9. 

A comment on the Lamy and Perrin volume scattering function is also in order. Lamy and 
Perrin (1986) adopt a volume scattering function (which is proportional to the product aa\ip 
used here) that varies with heliocentric distance as r~ Va where v c is chosen so that the dust cross 
sectional density cr(r) falls off as r~ Ul where V\ = 1. Although Lamy and Perrin (1986) provide 
excellent arguments to motivate their approach, their assumption will not be adopted here, but 
only because we wish to compare our findings to other studies of the zodiacal cloud that similarly 
do not apply this assumption. However our findings reported below are easily recalibrated if the 
Lamy and Perrin volume scattering function is preferred. If it is assumed that does indeed vary 
as r~ Uc then this implies that the dust albedo also varies as the power law a(r) = a\(r jr\)~ Vc where 
a\ is the dust albedo at r = 1 AU. Accounting for this alternate interpretation thus requires the 
substitutions a — > a± and v — > v\ + u c in Eqs. (3-4), but this has no substantive effect on their form. 
The only significant changes to our findings would be (i.) that the albedo quoted in Section 4.3 
should be interpreted as the albedo a\ for dust at 1 AU and that (ii.) the integrated dust surface 
densities and masses given in Eqs. (10-13) need to be reevaluated for the alternate dust density 
power-law v\. But as long as v c (which might be as large as 0.45) is smaller than unity then these 
revisions will change our finding by factors that are also of of order unity. 



4.2. Vertical variations 

Estimates of the dust grains' vertical distribution is obtained by first developing a simple yet 
plausible model for the dust inclination distribution g{i). The dust latitude distribution h([3) is 
then calculated from this inclination distribution, and a surface brightness map of the model dust 
cloud is generated using the line-of-sight integral Eq. (3). Then by comparing isophotes of the 
model cloud to the observations (Fig. 5), the suite of dust models that are consistent with the 
observations are readily obtained below. 
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4-2.1. inclination distributions 

Section 4.2.2 will describe the various inclination distributions Qj{i) for distinct dust-sources 
i.e., the asteroids, comets, and interstellar sources that are indicated by the j subscript. The 
relationship between population j's latitude distribution hj(/3) and its inclination distribution gj(i) 
is 

9j{i)di 



= / /■ ■ 2 R W 

J/3 v sir i - sur 

(Divari 1968; Divine 1993; Brown 2001). Although Eq. (6) is formally derived for bodies in circular 
orbits, it nonetheless provides reliable results even for bodies in very eccentric orbits [c.f. Brown 
(2001)]. Note that an isotropic cloud has an inclination distribution gj(i) = (2/tt) sin i which results 
in a latitude distribution hj((3) = 1. Accordingly, the total distribution will be written as a sum 
over several possible components, one being an isotropic source with an iso subscript and the N 
other sources having a gaussian distribution of inclinations: 

... 2 . . [ 1 j = iso 

yjw vr | Cj e- {l l u ^ I 2 otherwise v ' 

where <7j is the standard deviation of each component's inclination distribution and cj is a normal- 
ization constant such that hj (0) is unity. If each population's cross-sectional density is assumed to 
vary as a distinct power-law oc r~ v \ then the zodiacal light's total surface brightness becomes [see 
Eqn. (3)] 

z(M) = E/^-(^) = aa ^ (—) B oE^r r^)M/%))si^%)^- (8) 

where the coefficients fj indicate each population's relative contribution to the cross-sectional 
density in the ecliptic at r\ = 1 AU. Of course each population may also have a distinct albedo 
aj, in which case the fj in Eq. (8) should be replaced by (aj/a)fj where a is the 'effective' albedo 
that would be obtained if only a single dust population was responsible for all of the zodiacal light. 
However the individual albedos of asteroidal and cometary dust are not known with any certainty 
so we will simply set aj/a = 1. 



4-2.2. source populations 

The inclination distribution for the meteoritic dust complex is related to the inclinations of 
the sources of that dust, namely, colliding asteroids, active dust-producing comets, and interstellar 
dust sources. The upper portion of Fig. 10 shows the inclination distribution of near-Earth, main- 
belt, and Trojan asteroids. The smooth curve in this figure also shows that the observed inclination 
distribution can be qualitatively represented by the form g(i) oc sin ie~ 

-(i/a) 2 /2 with a = g^o^ Not 

however, that this curve seriously undercounts asteroids with inclinations i > 20° that are about 
6% of the total asteroid population. 



-14- 



The inclination distributions for comets having perihelia q < 2.5 AU that are presumably active 
producers of dust are also shown in Fig. 10 where they are divided into three dynamical classes: 
the Jupiter-family comets (JFCs), the Halley-type comets (HTCs), and the Oort Cloud comets 
(OCCs). The JFCs have the lowest inclinations since they likely originated in the (relatively) 
low-inclination Kuiper Belt (Levison and Duncan 1997). Conversely, the wide-ranging OCCs have 
the highest inclinations since their orbital planes have been randomized by the galactic tide and 
passing stars (Duncan et al. 1987). However the HTCs have intermediate inclinations as they likely 
originate in a somewhat flatted inner Oort Cloud (Levison et al. 2001). The smooth curves in Fig. 10 
also show representative inclination distributions with a = 8° for the JFCs, a = 33° for the HTCs, 
and g(i) oc sini for the isotropic OCCs. Interstellar dust should also be distributed isotropically in 
the inner solar system and thus have a similar inclination distribution. Note, however, that these 
cometary ex's should be regarded as rough estimates since the apparent orbital distributions suffer 
from selection effects that tend to favor the discovery of comets in low-perihelia, low-inclination 
orbits. In fact, the underabundance of high-inclination OCCs in Figure 10 suggests that selection 
effects may be especially severe for these single-apparition comets. 

Evidently, there are three classes of dust sources having distinct inclination distributions: a 
lower inclination population composed of asteroids and JFCs having a ai ow ~ 7° , a higher inclination 
population composed of HTCs having a (Thigh — 33°, and an isotropic population of dust from OCCs 
and interstellar sources. With these <j's in hand, the normalization coefficients Cj appearing in Eq. 
(7) are obtained by numerically integrating Eq. (6) and requiring hj((3 = 0) equal unity, which 
yields ci ow = 10.27 and c high = 2.190. 

The remainder of this study shall assume that the dust observed in Fig. 5 have the same 
inclination distributions as their source populations. However this need not necessarily be true since 
dust grains will slowly spiral sunward due to Poynting-Robertson (PR) drag. As a consequence of 
this radial mobilization, dust grains can traverse secular resonances with the planets where they 
can experience additional inclination excitation. The degree of this excitation is size dependent 
since a larger particle will drift at a slower rate across a resonance and thus experience greater 
inclination excitation. Dynamical models indicate that grains with radii R > 100 /urn will suffer 
significant inclination-pumping, i.e. Ai > 10°, whereas grains smaller than R < 30 fim suffer 
only modest pumping, i.e. Ai < 3° (Jackson and Zook 1992; Dermott et al. 2001; Grogan et al. 
2001). However the cross-sectional area of interplanetary dust is dominated by grains having radii 
10 ^m < R < 100 /urn which peaks at R ~ 30/t/m (Grim et al. 1985), so this additional source of 
inclination-excitation is at most a marginally important effect that is not modeled here. But even 
when a dust grain is far from a secular resonance, the planets secular perturbation will still excite 
forced inclinations that are of order a few degrees. These latter perturbations, which are also not 
treated by our model, introduces a warp in the zodiacal cloud whose inclination varies with the 
dust semimajor axis. 



-15- 



4.3. The abundance of asteroidal and cometary dust 

With the inclination distributions characterized, Eq. 7 can be inserted into Eq. 6 and the 
latitude distributions hj(/3) can be evaluated numerically for each population. This in turn allows us 
to compute maps of each population's surface brightness Zj(0, 4>) on a geocentric latitude-longitude 
grid using Eq. 8, and isophotes for example populations are shown in Fig. 11. A synthetic surface 
brightness map of the zodiacal light is formed by selecting the power-laws Vj for each population, 
computing the Zj(6,(p) maps, and then coadding the maps with proportions fi ow , fhi ,and fi so . 
Although the synthetic map has seven parameters, i.e, three fj, three Vj, and ao\r\, the available 
parameter space is quite limited since < f 3 < 1. Also recall the single power-law evident in 
the ecliptic surface brightness profile, Fig. 8. The absence of a broken power-law in the ecliptic 
profiles suggests either (a) all populations have a similar Uj ~ 1.45 radial variation, or (b) there 
is a single population having Uj ~ 1.45 that contributes most of the observed light. Note that the 
asteroidal dust likely varies as r~ l which is expected for dust that spirals into the field of view via 
PR drag 2 . However the cometary components should vary faster than r~ l since these dust grains 
can are produced in situ at rates that also vary with distance r. Consequently we expect the three 
populations to have 1 < Vj < 2 or so. 

We have scanned the Vj parameter space between Uj = 1.0 to Uj = 2.5 in increments of 
about Auj = 0.5. It is straightforward to scan the remaining fj parameter space for a given Vj 
triplet, and in general there is only a single set of fj parameters that agrees with the observations 
in a least-squares sense. Best agreement with the data is achieved when the Vj and fj take 
the parameters listed in Table II with the error bars indicating the range of possible fits that are 
marginally acceptable. Contours for this model are shown in Fig. 12 which compares quite favorably 
to the observed isophotes. However a model having all the Vj = 1.45 and fi ow = 0.34 ± 0.06, 
fhigh = 0.51 ± 0.11, and fi so = 0.15 ± 0.06 yields isophotes that are very similar to Fig. 12 with 
agreement that is almost as good. This indicates that this model is particularly sensitive to the 
bright, high-inclination dust population that is distributed over a wide range of latitudes, but that 
the power-law variations in the low-latitude as well dimmer isotropic populations are less well 
constrained. In general, we find acceptable solutions only when 1.0 < vi ow < 1.45, fhigh — 1-45, 
and 1.45 < ui so < 2.0 having approximate abundances of fi ow ~ 0.4, fhigh ~ 0.5, and fi so ~ 0.1. 

Recall that there is still a population of high inclination asteroids having i as t ~ 22° that have 
not yet been considered (see Fig. 10); these asteroids represent about 6% of the total asteroid 
population. Might these asteroids be a significant source of the high latitude dust seen in Fig. 
5? If so, then their fractional contribution to the ecliptic surface brightness would be of order 
fast ~ 6%fiow (J iow/'i'ast ~ 0.01 when smeared out over an annulus that is i a st/&iow ~ 3 times thicker 
than that inhabited by the lower-i asteroids. Consequently, these high-i asteroids are not likely to 



2 This assertion is true when the dominant light-reflecting grains have an orbital lifetime due to PR drag that is 
shorter than their collision lifetimes, which is the case for grains having radii smaller than about ~ 100 /im (Griin et 
al. 1985). 
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be a significant source of high-latitude dust since their ecliptic contribution is only ~ 2% that of 
the high-inclination population's contribution fhigh- 

The inferred dust latitude distribution h((3) is also shown in Fig. 13 as well as the weighted 
contributions fjhj(/3) by the low, high, and isotropic dust populations. This figure shows that at 
ecliptic latitudes (3 > 15°, more than 90% of the cross-section is contributed by dust that are in 
comet-like orbits (e.g., HTCs and OCCs). The spatial distribution of the dust cross section is 
also shown in Figure 14. Although these contours are rather similar to that inferred by Kelsall et 
al. (1998) from the COBE observations, they do differ in detail due to the different assumptions 
built into each model. Specifically, the Kelsall et al. model employs an empirical function (namely, 
a modification of the familiar fan model) to describe the dust latitude distribution h{0) of their 
'smooth cloud' (which is the principle component of that model), whereas we allow for three distinct 
distributions hj((3) that are instead based upon known comet and asteroid inclination distributions. 
Although both modeling efforts adopt very different treatments of the dust latitude distributions 
h(f3), and these models were also applied to data acquired at rather different wavelengths (optical 
versus near and far infrared), the inferred dust density distributions reported here and by Kelsall 
et al. (1998) are quite similar. The density distribution given in Fig. 14 is also reminiscent of the 
familiar fan model that assumes a(r,(3) = a\{r /r\)~ u e~ k \ sin ^l. However a parameterization of this 
form provides at best only a qualitatively correct estimation of the density map of Fig. 14 when 
k ~ 1.5; it still fails to reproduce this Figure in detail, especially at high latitudes beyond (3 > 45° 
and at distances beyond r > 2 AU. 

The inferred inclination distributions for all of the dust populations are also shown in Fig. 15. 
We also note also that these distributions are rather similar to that reported by Divine (1993) who 
inferred dust size and orbital distributions from a wide suite of dust observations (e.g., microcraters 
on lunar samples, spacecraft dust-impact experiments, as well as other zodiacal light observations). 
However Divine did not comment on the implications of this inclination distribution, which we 
regard as one of the more interesting findings of this study. 

Since &j(r, f3) = j jcr\(r / l r\)~ v i hj(/3) is the density of population j's dust cross-section, its total 
dust cross-section contained within a sphere of radius r is obtained from the volume integral 



where the latitude integration is 7j = J* hj(f3) cos f3d/3, which numerically integrates to 7; D?i , = 
0.143, jhigh = 0.619, and 7j so = 1.00. Adopting the best-fitting Vj and fj parameters given in Table 
II, each population's total dust cross-section interior to r± = 1 AU is Sj(ri) = {0.032, 0.200, 0.050}x 
4-Kairf, respectively, for the low, high, and isotropic populations. Although the low-inclination dust 
from asteroids and JFCs contributes fi ow = 45% of the dust cross-sectional density in the ecliptic, 
at least 89% of the dust interior to a r\ = 1 AU sphere is contributed by sources in comet-like 




(9b) 



(9a) 
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orbits (e.g., HTCs and OCCs plus an unknown fraction from JFCs). 

Using data from spacecraft dust collection experiments as well as studies of lunar microcraters, 
Grim et al. (1985) estimate the spatial density of the ecliptic dust cross section at 1 AU to be 
<7i ~ 4.6 x 1CT 21 cm 2 /cm 3 . Since Table II reports ao\r\ = (7.4 ± 1.8) x 1CT 9 , this implies that the 
dust have an effective albedo of a ~ 0.1. Note, however, that these estimates for o\ and thus a are 
probably uncertain by a factor of ~ 2 since these dust collection experiments largely measure dust 
fluxes versus particle energy, and that their conversion to a dust cross-section requires assumptions 
about the dust velocities and their bulk densities. 

Summing Eq. (9b) over all populations yields the total dust cross-section contained within a 
sphere of radius r: 



where each term gives the contribution by the low i, high i, and isotropic populations, respectively. 
Note that this expression only applies interior to the dust-producing portion of the asteroid belt, 
i.e. interior to r ~ 3.3 AU (Hanner et al. 1974). If we consider a sphere of radius ri = 2 AU 
enclosing the orbits of the terrestrial planets, then S(r2) = 1.6 x 10 10 km 2 , which is about 50 
times the total cross-section of the terrestrial planets. This estimate illustrates one of the main 
difficulties challenging efforts to detect extra-solar planets via direct imaging or interferometry at 
optical wavelengths: if terrestrial extra-solar planets are also embedded in solar system-like dust, 
then one will need to resolve planetary systems to fairly small spatial scales in order to discriminate 
the starlight reflected by planets from that reflected by dust. 

The Pioneer 10 spacecraft detected asteroidal dust out to a heliocentric distance of r% ~ 3.3 
AU (Hanner et al. 1974), so Eq. (9b) indicates that the low-inclination dust component has a total 
surface area £; ou ,(r3) = 6.8 x 10 9 km 2 . Note that observations of the IRAS dust bands serve as a 
comforting reality check on our findings since the total surface area associated with the three most 
prominent asteroidal dust bands is £& a n<2 = 4.7 x 10 9 km 2 (Grogan et al. 2001); the remaining dust 
must then be due to other minor asteroid families, non-family asteroids, and JFC comets. Also, 
if the light-reflecting dust seen in Fig. 5 can be attributed to grains having a characteristic radius 
R c , then the total number of grains interior to distance r is Nj(r) ~ £j(r)/7T-R 2 and their enclosed 
mass Mj(r) ~ ApR c T,j(r)/3 is 



where p is the grains' bulk density In this case the total mass of the light-reflecting component of 
asteroidal dust is at most of order 




(10) 




(11) 




(12) 
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The interplanetary dust mass distribution peaks at R c ~ 100 /j,ia (Griin et al. 1985), so the above 
dust mass-limit is equivalent to an asteroid that is about 12 km across. Of course this limit is 
valid only if our model, which is based on observations of dust orbiting at r < 0.6 AU, can be 
reliably extrapolated out to the asteroid belt. Nonetheless, the similarity between our dust model 
and that inferred from the COBE observations, which are sensitive to dust in the 1 < r < 3 AU 
interval (Kelsall et al. 1998), indicate that our extrapolation is indeed valid. We also note that the 
mass limit obtained here about 3.5 times the mass of the asteroidal dust bands detected by COBE 
(Reach et al. 1997), but keep in mind that our limit is also contaminated by dust from JFCs. 

The isotropic cloud of dust is also quite interesting, and the following discussion assesses the 
relative dust contribution from Oort Cloud comets versus interstellar sources. The dust seen in 
Fig. 5 having a typical elongation of e ~ 15° orbit at a heliocentric distance of r ~ r\ sine ~ 0.3 
AU, and this dust has a cross-sectional density of o"j so (0.3 AU) = /j so <Ti(0.3 AU/ri)~^ so ~ 3 x 
10 -21 cm 2 /cm 3 in the ecliptic. The interstellar fraction is inferred from the flux of interstellar 
dust measured by impact detectors onboard the Galileo and Ulysses spacecraft. These detectors 
measured an interstellar dust flux of /* = 1.5 x 10 -8 grains/cm 2 /sec [after correcting a typo in 
Griin et al. (1997)]. These interstellar grains have a mean mass of m+ ~ 3 x 10~ 13 gm, so their 
characteristic radius is ~ 0.3 fim. It will be assumed here that the interstellar dust flux is roughly 
constant throughout the solar system since radiation pressure roughly balances solar gravity for 
grains of this size. Since interstellar matter approaches the solar system with a velocity-at-infinity 
of Voo = 25 km/sec (Frisch 2000), the number density of interstellar dust is n* ~ /*/i>oo ~ 6 x 10~ 15 
cm -3 and their cross-sectional density is a+ ~ 7ri? 2 n + ~2x 10~ 23 cm 2 /cm 3 . This is only about 1% 
of the observed cross-sectional density (Ji so , which indicates that the isotropic portion of the dust 
seen in Fig. 5 comes predominantly from Oort Cloud comets. 

Oort Cloud comets have semimajor axes a ~ 10 4 AU, and those comets passing sufficiently 
close to the Sun will sublimate gas and dust that gets injected into orbits similar to their parent 
comets. If we naively extrapolate these wide-ranging dust grains out to Oort Cloud distances using 
the inferred Vi so = 2 radial power-law, the total Oort Cloud dust mass is 

1 * 1018 (life*) G£0 (ioTSu) < 13 > 

which has a equivalent to a ~ 30 km comet. However this mass is easily uncertain by orders of 
magnitude due to uncertainties in the size and bulk density of the dust as well as the radius of the 
Oort Cloud. Ultimately these distant dust grains will be stripped from the solar system as they are 
swept up by the interstellar gas and dust that flows through the solar system. Thus it is conceivable 
that the Sun also has a vast but tenuous tail of Oort Cloud dust. If so, this dust tail would be 
oriented in the downstream direction of the local interstellar flow which has a heliocentric ecliptic 
longitude, latitude of (74.7°, —4.6°) (Frisch 2000) or an equatorial right ascension, declination of 
(73.9°, +18.0°). 
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5. Summary and Conclusions 

Using the Moon to occult the Sun, the Clementine spacecraft used its navigation cameras 
to map the inner zodiacal light at optical wavelengths over elongations of 3 < e < 30° from 
the Sun. Since the zodiacal light is sunlight that is reflected by interplanetary dust, this map 
provides a measure of the dust grains' radial and vertical variations spanning heliocentric distances 
of 0.05 < r < 0.6 AU, i.e., from about 10 solar radii to just interior to Venus' orbit. The integrated 
zodiacal light seen over the 60° x 60° field of view has a visible magnitude my = —8.5, indicating 
that the meteoritic complex is one of the brightest members of the planetary system, second only 
to the full Moon. 

The averaged ecliptic surface brightness of the zodiacal light falls off as Z(e) oc e ~ 2 - 4 5±o.05 
which suggests that the dust cross-sectional density nominally varies as a(r) oc r - 1A 5±o.05 
see below). This surface brightness also indicates that the dust obey ao\r\ = (7.4 ± 1.8) x 10~ 9 . 
Assuming that the dust have an ecliptic cross-sectional density of o\ = 4.6 x 10~ 21 cm 2 /cm 3 at 
r± = 1 AU (Griin et al. 1985), this implies that the dust have an effective albedo of a = 0.1 that is 
perhaps uncertain by a factor of 2. Asymmetries of ~ 10% in the zodiacal light's surface brightness 
are evident at elongations e ~ 15° in directions east-west as well as north-south of the Sun, and 
these asymmetries may be due to the giant planets' secular gravitational perturbations. However 
simple digital filtering of the data (i.e., unsharp masking) does not reveal any other subtle features 
such as dust bands associated with asteroid families or dust trails associated with individual comets. 

In order to assess the relative contributions to the interplanetary dust complex by asteroids 
and comets, we have modeled the zodiacal cloud as being due to three dust populations having 
distinct inclination distributions. One dust population is assumed to have low inclinations that 
are distributed as a gaussian with a standard deviation oi ow — 7° that is characteristic of both 
asteroids as well as Jupiter-Family comets (JFCs). A higher-inclination population corresponding 
to the Halley-type comets (HTCs) is assumed to have a (Thigh — 33°, and the third population is an 
isotropic cloud of dust from Oort Cloud comets. This simple model is applied to the observations 
and very good agreement is achieved for the parameters listed in Table II (see Fig. 12). It should 
be noted that the best-fitting model has a radial power-law v\ ow = 1.0 for the asteroidal + JFC 
population, which is consistent with dust delivery via Poynting-Robertson (PR) drag yet shallower 
than the nominal v = 1.45 power-law quoted above. Also, the dust from HTC appears to follow 
a Vhigh = 1-45 power-law while the Oort Cloud dust varies as z/j so = 2.0. Interstellar dust also 
contributes to this isotropic cloud, but only at the ~ 1% level. Yet despite this mixture of power- 
laws, the resulting surface brightness profile still varies close to the observed Z(e) oc e~ 2 ' 45 . Note, 
however, that acceptable agreement with the data is also achieved when all populations have 
Vj = 1.45; see Section 4.3 for the allowed range of model parameters. 

The best-fitting model indicates that about fi ow = 45% of the dust cross-section in the ecliptic 
at r\ = 1 AU comes from asteroids and JFCs. But when a 1 AU-radius sphere is considered, at 
least 89% of the integrated dust cross section comes from sources in comet-like orbits. However 
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it should be noted that these findings are inferred from a 'static' model of stationary dust grains. 
This rather simple approach is applicable provided dust grains (and in particular, asteroidal dust) 
do not experience substantial-inclination pumping as they evolve sunwards due to PR drag and 
cross orbital resonances with the planets. Although this appears to be the case for grains smaller 
than ~ 100 fim that are the dominant source of reflected sunlight, inclination pumping is certainly 
of greater importance for the larger dust grains that can drift across resonances at slower rates. 
In this instance, 'dynamic' models that include radiation and gravitational forces are preferred 
as they can faithfully follow a dust grain's orbital evolution from source to sink, and these more 
sophisticated models have been used to extract dust properties from the infrared observations of the 
outer zodiacal light acquired by IRAS (cf. Dermott et al. (2001); Grogan et al. (2001)). Of course 
these dynamic dust models are also applicable to the optical observations of the inner zodiacal 
light examined here, and we suspect that any effort to simultaneously fit a dynamic dust model to 
the full suite optical and infrared IRAS, COBE, and Clementine observations should yield an even 
more tightly constrained picture of the interplanetary dust complex spanning a very wide range of 
heliocentric distances 0.05 < r < 3.3 AU. And in order to facilitate any such effort, the Clementine 
map of the inner zodiacal light is available from the authors by request. 

The inferred cross-section of dust orbiting interior to 2 AU is about 50 times that of the 
terrestrial planets. This suggests that any effort to directly detect terrestrial extrasolar planets at 
optical wavelengths is faced with the daunting task of distinguishing the faint starlight reflected 
by such planets from the far brighter signal anticipated from any exozodiacal dust. When these 
results are extrapolated out into the asteroid belt, the total mass of the light-reflecting asteroidal 
dust component is at most ~ 2.3 x 10 18 gm. Note that this mass-limit is also contaminated by 
dust from JFCs, and it corresponds to an asteroid that is at most ~ 12 km across. 

When these results are extrapolated out to Oort Cloud distances of a ~ 10 4 AU, the inferred 
mass of Oort Cloud dust is ~ 10 19 gm (but uncertain by orders of magnitude), which is equivalent 
to a 30 km comet. This dust is ultimately stripped from the Sun by the interstellar gas and dust 
that flows around and through the solar system. This then suggests that the Sun, and perhaps 
also other stars having cometary Oort Clouds, each have vast but tenuous stellar dust tails that 
are oriented in the downstream direction of the local interstellar flows. 



- 21 - 



A. Appendix A 



To relate a pixel's (x, y) coordinates to the right ascension and declination (a, 5) it subtends 
on the sky, place a Cartesian coordinate system with its origin at the center of the lens with the x, 
y axes parallel to the CCD's rows and columns (see Fig. 1). The z axis is the camera's optical axis 
which intercepts the CCD at the pixel having coordinates (x ,y ) = (191,286). Let the vector r* 
point to a star having angular coordinates (6,<p); these will be called the lens coordinates. In carte- 
sian (x,y,z) components, the star's lens coordinates are = (R sin (ft cos 9, R sin 9, R cos (ft cos 9), 
and the length of this vector is chosen to be the lens radius R. An image of this star also forms at 
the back side of the lens at — r*. Fig. 1 shows that the fiber optic pipes this starlight to a spot on 
the CCD that lies a distance X away from the optical axis in the x direction and Y away in the 
— y direction. Note that the orientation of the X and Y axes are chosen to point in the customary 
manner such that when y points to equatorial north, the +Y direction is north and the +X direc- 
tion is west. If I is the physical size of a pixel, then X = l(x — x Q ) = — r+ ■ x, Y = l(y — y Q ) = r* ■ y, 
and 



relates the star's (x,y) coordinates on the CCD to its lens coordinates (9, (ft). The plate-scale 
p = i/R is simply the ratio of the pixel width to the lens radius, and is also the angle subtended 
by the pixel at the optical axis. 

Now relate the lens coordinates (9, (ft) to equatorial coordinates (a, 5). Suppose the optical 
axis z points to right ascension a Q and declination 5 , and the y axis differs from equatorial north 
by a rotation about the z axis by angle r. This angle is the position angle of equatorial north, and 
on the CCD it is measured from the Y axis towards the east. In this coordinate system the star's 
cartesian components are r*(a, S) = (R sin a cos 5, R sin 5, R cos a cos 5) when expressed in terms of 
its equatorial coordinates. It is straightforward to show that this star's position vector r^(a, S) in 
equatorial coordinates can be obtained from the its position vector in lens coordinates r*(</>, 9) after 
performing the following rotations upon the lens coordinate system: 



where lZi(u) is the matrix that rotates a right-handed coordinate system about axis i by angle uj. 
Performing the rotations and equating the (x, y, z) components in Eq. (A2) yields the relationship 
between a pixel's lens coordinates (9,(j)) and its equatorial coordinates (a, 5): 



x — x = — sin (j) cos 9 /p 
V~Vo = sin9/p 



(Ala) 
(Alb) 



r*(a,<5) = K y {-a )n x {5 )K z {-T)T±{9,4>) 



(A2) 



tan a 



sin 5 



g cos 5 + h sin S 
k sin a Q + f cos a 



k cos a a — f sin a Q 



(A3a) 
(A3b) 



where f,g,h, and k are shorthand for 



/ = cos r sin <j) cos 9 — sin r sin 9 



(A4a) 
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g = sin r sin <f> cos 9 + cos r sin 9 (A4b) 
h = cos(/>cos# (A4c) 
= hcos5 Q — gsin5 D . (A4d) 



Next, get the star's geocentric ecliptic longitude and latitude (A, 0) where 
r*(A, 0) = (i?sin A cos 0, i?sin 0, Rcos A cos 0). These are obtained by rotating the ecliptic co- 
ordinates about the z axis by the Earth's obliquity a = 23.4393°, i.e., r*(A, 0) = lZ z (a)r*(a, 8). 
Equating cartesian components yields 

sin = cos a sin 5 — sin a sin a cos 5 (A5a) 

cos a sin a cos 8 + sin a sin 5 , . „, . 

tan A = . A5b 

cos a cos o 

This gives the star's longitude measured from the vernal equinox. To get the star's longitude 
relative to the Sun, use Eq. (A5b) to compute the Sun's longitude A Q and form the difference 
A-A . 

The observed angular separation $oi between two known stars can be used to determine the 
detector's plate-scale p. If the stars have lens coordinates (9o,(po) and (9i,(/>i), which depend on 
p, and equatorial coordinates (ao,5o) and (ai,6i), then spherical geometry gives their angular 
separation: 

cos $oi = cos #o cos 9\ cos(0o — 4>i ) + sin #o sin 9\ (A6a) 
= cos 8q cos 8i cos(ao — a\) + sin8o sinJi. (A6b) 

This equation can be solved for p when coupled with Eqs. (Al). Another useful quantity is the 
angular separation between star and the optical axis at (9i,(/)i) = (0,0), 

cos <I>o = cos #o cos 4>q. (A7) 

Using Fig. 1, it can also be shown that an off-axis pixel subtends a larger solid angle f2(<I>o) than 
the on-axis pixel: 

- 1 

flo cos$ 

where = P 2 is the solid angle of the pixel at the optical axis. For example, a pixel in the CCD 
corner at <E>o = 25° sees a solid angle that is 10% larger than the on-axis pixel. 



B. Appendix B 



Flatfield data for camera B were acquired prior to launch in a laboratory at Research Support 
Instruments, which is the subcontractor that integrated the star trackers into the Clementine 
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payload. These images were acquired with star tracker B exposed to an integrating sphere, which 
is a spherical light source having a uniform radiance. The averaged image acquired in the lab is 
shown in Fig. 16A. The following describes several problems in these data, as well as the procedures 
used to construct the star tracker flatfield from these data. 

(i.) The integrating sphere was too bright to produce unsaturated images at the star tracker's 
shortest possible exposure time. To get unsaturated images, a Wrattan neutral density ND2 filter 
was placed in front of the camera. It is thus conceivable that if there are any spatial variations 
8T(x,y) in this filter's transmission coefficient T(x,y), the resulting flatfielded data would then 
exhibit fractional errors 5T/T. However we have ruled out this possibility by examining distinct yet 
overlapping zodiacal light fields acquired during different orbits with different camera orientations. 
In every case, the photometry in each overlapping, flatfielded image was self-consistent. This 
indicates that any errors 5T(x, y) that may have been introduced into the flatfield by the use of 
the neutral density filter are negligible. 

(ii.) The dark current was not measured when the flatfield data were acquired. This could 
have been quite problematic since the the subtraction of the dark current has a dramatic affect 
upon the brightness of the flatfield center relative to its edges (see Fig. 16). However a good flatfield 
is still recoverable due to the fortunate presence of the five round blemishes seen in Fig. 16. These 
blemishes correspond to a ~ 10% decrease in the camera's sensitivity, and their cause is unknown. 
These blemishes are present in all of the raw zodiacal light images acquired by star tracker B 
(c.f. the lower right corner of Fig. 4A). However the relative depth of the blemish increases with 
increasing dark current d, so it is straightforward to find the appropriate value for d that yields a 
flatfield that removes the blemishes from all of the zodiacal light images. Note that the flatfielded 
image in Fig. 4B does not exhibit these blemishes. 

(to.) Vignetting by the mount that secured the neutral density filter to the camera also 
darkened the outer edges of the flatfield data (see Fig. 16A). Interior to the vignetted region, the 
flatfield's surface brightness empirically varies as a simple polynomial 

flat(x,y) = l-1.95$g(a;,y) (Bl) 

that is a function of the angular distance <I>o of pixel (x,y) from the the optical axis (see Eq. A7). 
In order to reconstruct the flatfield at the edges, Eq. (Bl) is extrapolated into the vignetted region 
at $o > 20°, which results in the flatfield shown in Fig. 16B. Gaussian noise is also added to the 
extrapolated region in amounts comparable to that seen just interior. 

The pixel-to-pixel variations seen over short spatial scales in Fig. 16 is simply photon-counting 
noise — they do not represent real changes in the camera's sensitivity. This becomes evident when 
flatfielding these images. Ordinarily, flatfielding an image makes it look 'cleaner' since the pixel- 
to-pixel variations in the CCD's sensitivity are removed from the image. However the application 
of the flatfield seen in Fig. 16B actually makes the zodiacal light images look a little bit noisier. It 
may be concluded that the short-wavelength variations seen in the flatfield are simply noise, and 
that we are justified in adding gaussian noise to the extrapolated parts of the flatfield. (Failure to 
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add this noise creates a lower-noise zone at the edge of a flatfielded image having a very artificial 
appearance.) The standard deviation of this noise near the optical axis is about 2% of the flatfield's 
surface brightness, and it increases to about 7% at the edge. 

Figure 16 and Eq. (Bl) show that the surface brightness of the flatfield decreases with distance 
from the optical axis by as much as 40% at the CCD corners where <J>o = 25°. This seemingly runs 
contrary to expectations since Appendix A shows that pixels further from the optical axis subtend 
a solid angle that is larger by a factor of 1/ cos(<I>o) (Eq. A8), and perhaps should be brighter. 
However this effect is more than offset by transmission inefficiencies in the fiber-optic and also 
vignetting by the lens housing, both of which get more severe for light entering the lens at larger 
angles (Lewis et al. 1991). Consequently, dividing a raw image by this flatfield simultaneously 
corrects for two effects: it compensates for the camera's reduced detection efficiency at higher <l?o 
as well as for the pixels' larger solid angle. 

Because the reconstructed flatfield of Fig. 16B compensates for a pixel's larger solid at higher 
$0; it should only be used when doing photometry on distributed sources such as the zodiacal light. 
But if photometry on point-sources is desired, a "point-source" flatfield, which is just Fig. 16B 
x cos(<l>o), must instead be used to flatfield the images. 

The observed intensity of several stars imaged during orbit 66 are used to test this reconstructed 
flatfield. The spacecraft was maneuvering while this image sequence was acquired, so stars seen in 
different images wander across the field. Figure 17 plots the normalized intensity of five bright stars 
seen in different images versus their angular distance $o from the optical axis. The grey curves 
are their intensities prior to flatfielding, which shows how the camera's detection efficiency drops 
with $o- The black curves shows the intensities of these stars when the point-source flatfield is 
used. This test demonstrates that when the reconstructed flatfield is used, the stars' intensities arc 
constant when imaged at different position on the CCD, as they should. 

This paper is dedicated to the memory of Herb Zook, without whom none of this would have 
been possible. The authors also thank Paul Spudis who, as the Clementine deputy scientist, allowed 
the spacecraft's navigation cameras to be pointed sunwards. The authors also thank Julie Moses 
for helpful discussions on interactions between interplanetary dust and the interstellar medium and 
Ronna Hurd for composing Fig. 1. The authors also thank William Reach and an anonymous 
reviewer for comments that led to several improvements in this paper. This paper is contribution 
#1117 from the Lunar and Planetary Institute which is operated by the Universities Space Research 
Association under NASA contract NASW-4574. 
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Table I. Zodiacal Light Observations 



Orbit 


Date 


Ecliptic 


Total 


number 


(1994) 


longitude 


exp. time 




(degrees) 


(seconds) 


66 


March 5.9 


345.1 


2.6 


110 


March 15.0 


354.2 


0.3 


110 


March 15.0 


354.3 


0.6 


164 


March 26.4 


5.5 


8.7 


193 


April 1.5 


11.5 


11.6 


206 


April 4.0 


14.0 


11.4 


253 


April 13.9 


23.7 


13.9 



Table II. Parameters for the Best Fit 



3 


v 3 


n 


low 


1.00 


0.45 ±0.13 


high 


1.45 


0.50 ±0.02 


iso 


2.00 


0.05 ±0.02 




cunn = (7.4 ± 1.8 


) x 10- y 




Fig. 1. — A schematic of the star tracker camera which has three principle components: a spherical 
lens and a fiber optic that pipes light from the backside of the lens to the CCD detector. 
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4000 6000 8000 10000 

wavelength (A) 

Fig. 2. — The dots indicate the relative sensitivity of the star tracker camera to light of discrete 
wavelengths of A = 5000 to 9000 A sampled over 1000 A intervals; these data are provided by J. F. 
Kordas (private communication). The dark curve simply connects the dots and also assumes that 
the camera sensitivity is zero at A = 4000 and A = 10000 A. The grey box is an equivalent square 
bandpass. 
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ENP 

GNP 




Fig. 3. — The narrow lines show the heliocentric ecliptic longitudes of the spacecraft's lines-of-sight 
through the inner zodiacal light, a few of which are labeled by their orbit number (see Table I). 
Longitudes are measured counter-clockwise from the direction of the vernal equinox 7. The dark 
ellipse is Earth's orbit about the Sun with ENP indicating the ecliptic north pole. The dashed 
circle represents the i = 3° tilt of the zodiacal light's midplane which has a longitude of ascending 
node 12 = 87° (Leinert et al. 1980). The galactic north pole GNP is also indicated, and it has an 
heliocentric ecliptic longitude X g = 180.0° and latitude [3 g = 29.8°. 
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Fig. 4. — Star tracker images shown at different stages of analysis. A. A typical raw star tracker 
image acquired during a 0.4 sec exposure in orbit 193. We saturate this linear greyscale at a surface 
brightness of 4 x 1CT 12 B in order to reveal the faint artifacts present in the data. The CCD's 
columns run left -right and ecliptic north/east are approximately up/left. The bright object left of 
the Moon is a saturated Venus whose signal has bled into the adjacent pixels. B. This logarithmic 
greyscale shows the master image for orbit 193 in the surface brightness interval 4 x 1CP 13 < Z < 
8 x 10 -11 B©. C. This cosmetically enhanced master has had its data-gaps filled in with images 
of the Moon and Venus. 




Fig. 5. — A mosaic of seven fields of the inner zodiacal light observed by the Clementine star tracker 
camera. The colorbar indicates surface brightness in units of 1O _13 B0. Ecliptic north and east 
are up and left in this mercator projection, and the field of view is 60 x 60°. Black indicates gaps 
in the data, and the Sun is drawn to scale at the center of the mosaic. Regions beyond ~ 10° 
northwest of the Sun are polluted by scattered light, and the "dimple" 20° east of the Sun is a 
lower signal/noise patch that was polluted by Venus. 
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Fig. 6. — The viewing geometry for light scattered by a dust volume element a distance r from the 
Sun at heliocentric latitude (5 measured from the ecliptic x-y plane. The observer has a heliocentric 
distance r\ and the line-of-sight (LOS) to the volume element a distance A away has a geocentric 
latitude 6, longitude (f>, and an elongation angle e measured from the sunward direction. The 
scattering angle <p is measured from the LOS to the anti-solar direction at the volume element. 
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scattering angle (p (degrees) 



Fig. 7. — The upper curve is the phase law obtained from the nominal volume scattering function 
of Lamy and Perrin (1986), while the lower curve is from Hong (1985). Note that the phase laws 
shown here are simply the volume scattering functions with the dust albedo a and cross sectional 
density a\ factored out and normalized to unity at ip = tt where the phase angle it — <p = 0. 
Although Hong (1985) gives a phase law for a v = 1 radial power law, we have used his Eqns. (10) 
and (14) to form the phase law shown here that is appropriate for a v = 1.45 dust distribution. 
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Fig. 8. — The surface brightness of the zodiacal light mosaic of Fig. 5 is plotted versus elongation 
angle e. Each profile is computed in a 10 degree-wide triangular aperture oriented North, South, 
East, or West of the Sun. The light-polluted field acquired during orbit 164 (which lies north of 
the Sun) is discarded before generating these profiles, as are pixels that subtend data-gaps, bright 
stars, or planets. 
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Fig. 9. — Ratios of the brightness profiles of Fig. 8 are plotted versus elongation angle e, where 
N/S and E/W indicate the north/south and east/west ratios. 
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Fig. 10. — The upper figure shows the inclination distribution N(i) for 70,383 asteroids having 
semimajor axes a < 6 AU brighter than absolute magnitude H = 15 (i.e, brighter than the 
completeness limit for asteroid surveys); this corresponds to asteroids having diameters between 
5 and 850 km. These data come from Edward Bowell's The Asteroid Orbital Elements Database 
obtained from the URL ftp://ftp.lowell.edu/pub/elgb/astorb.html. Also shown are the inclination 
distributions for 138 Jupiter-Family comets (JFCs), 27 Halley-Type comets (HTCs), and 223 Oort 
Cloud comets (OCCs). Only those comets having perihelia q < 2.5 AU listed in the Marsden and 
Williams (1999) catalog are used here. The smooth curves are g(i) oc sm(i)e~^^ a ^ I 2 with a = 6.2° 
for the asteroids, a = 8.0° for the JFCs, a = 33° for the HTCs, and g(i) oc sin(i) for the OCCs. 
Since this study of the zodiacal light is insensitive to prograde versus retrograde orbits, we have 
replaced each retrograde orbit having an inclination i > 90° with a prograde equivalent having an 
inclination 180° — i in order to improve the statistical significance of these figures. 
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Fig. 11. — Isophotes for the low- inclination dust population (e.g., dust from asteroids and JFCs), 
the high-inclination population (dust from HTCs), and the isotropic population (dust from OCCs 
and interstellar sources). These contours are computed using Eq. 3 and the Hong phase law shown 
in Fig. 7 and with the dust cross sectional density varying as a(r) oc r~ v i with vj = 1.45 for each 
population. The brightness of each isophote differs by a factor of 2. Isophotes for populations with 
Vj = 1 are shaped similarly but are shifted inwards while Vj = 2 isophotes are shifted outwards. 




Fig. 12. — The narrow black curves are isophotes of the zodiacal light map of Fig. 5, and the thick 
grey curves are isophotes for the best-fitting model having the parameters given in Table II. The 
brightest isophote corresponds to a surface brightness of 7.45 x 10~ n B and successive isophotes 
are dimmer by factors of 2. The outermost isophote is smoothed over a box 5 pixels or 0.38° wide. 
The rectilinear features seen in the outermost isophote are 'seams' in the mosaic image (see Fig. 
5); they are due to small errors in the offsets 5f that were removed from each image (see Section 
3). The jag in the innermost contour as well as the linear features seen at eft = —20° and —27° east 
of the Sun are all due to gaps in the data, and the bends in the isophotes north-northwest of the 
Sun are due to scattered light in the camera. 
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Fig. 13. — The total dust latitude distribution h((3) plotted versus heliocentric ecliptic latitude (5. 
Also shown are the weighted contributions by the low, high, and isotropic populations, fjhj(fi). 
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Fig. 14. — Contours of the inferred density of dust cross-section a(x, z) in cylindrical coordinates 
where x is the ecliptic distance from the Sun and z the height above the ecliptic plane. Adjacent 
contours indicate a factor of 1.5 change in the dust density and the dotted curve is where the dust 
density a(x,z) = o\ = 4.6 x 10~ 21 cm 2 /cm 3 . 
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Fig. 15. — The total inclination distribution g(i)/ sini plotted versus inclination i as well as the 
weighted contributions from the low, high, and isotropic populations. 
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Fig. 16. — A. The vignetted flatfield. B. The reconstructed flatfield. 
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Fig. 17. — The intensities of five stars (indicated by different plotting symbols) seen in different 
images acquired during orbit 66 are plotted versus each star's angular distance from the optical 
axis $o- All intensities are normalized to unity near <E>o = 0. The grey curves are extracted from 
images that are not flatfielded, and the dark curves are from images obtained with the reconstructed 
point-source flatfield. 



